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Abstract 

The normalized differential mean free path for volume scattering and the differential surface 
excitation probability for medium energy electrons travelling in Fe, Pd and Pt are extracted from 
Reflection Electron Energy Loss Spectra (REELS). This was achieved by means of a recently 
introduced procedure in which two REELS spectra taken under different experimental conditions 
are simultaneously deconvoluted. In this way, it is possible to obtain the unique reconstruction 
for the surface and volume single scattering loss distribution. The employed method is compared 
with a procedure that is frequently used for this purpose [Tougaard and Chorkendorff, Phys. Rev. 
B 35(1987)6570]. It is shown, both theoretically and through analysis of model spectra as well as 
experimental data that this method does not result in a single scattering loss distribution. Rather, 
it gives a mixture of surface, bulk and mixed scattering of any order. 
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1. INTRODUCTION 



Quantitative analysis of reflection electron energy loss spectra (REELS) is important for 
surface science both in a practical and fundamental sense since such spectra contain valuable 
information on the response of a solid to an external electromagnetic perturbation. From 
the fundamental point of view it is desirable to extract the frequency (ou) and momentum 
(q) dependent dielectric function e(q,u) from experimental data, being the characteristic 
quantity for the susceptibility of the solid to polarize under bombardment with electromag- 
netic radiation [1-5]. For practical applications of electron spectroscopy, information about 
the differential distribution of energy losses in an individual loss process is a prerequisite 
for quantitative analysis of electron spectra [6, 7]. Such data can be directly extracted from 
REELS spectra and converted to dielectric data by invoking linear response theory. 

The interaction of a reflected electron with a solid is quite complex since different types 
of collisions occur along the electron trajectory, and, furthermore, these scattering processes 
can occur repeatedly when the pathlength of the electron in the solid is long compared to 
the relevant mean free paths for scattering. The changes in the degrees of freedom of the 
electron are brought about by the interaction with the ionic subsystem, which mainly changes 
the direction of the particle in an elastic process and the interaction with the electronic 
subsystem, that gives rise to energy losses [7]. An energy loss process can be conceived as 
a decelaration of the probe electron by a polarization field that is set up in the solid by the 
response of the weakly-bound solid state electrons to the incoming electron [8] . This physical 
picture provides the connection between the dielectric function and the differential mean free 
path for inelastic electron scattering. Near a surface, the situation is further complicated 
through the occurence of so-called surface excitations [9-18]. These are additional modes of 
the inelastic interaction that result from the boundary conditions of Maxwell's equations near 
an interface between regions of different electromagnetic susceptibility. Surface excitations 
have a lower resonance frequency as compared to volume excitations and they occur on 
both sides of the solid-vacuum interface, but the probability for a surface excitation decays 
rapidly with the distance from the interface. 

The question then arises how to extract the single scattering loss distribution for volume 
and surface excitations from an experimental spectrum. An important step to resolve this 
problem was taken some time ago by the pioneering work of Tougaard and Chorkendorff 
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[6] who found a simple deconvolution formula for REELS spectra. However, the physical 
model on which it is based completely ignores surface excitations. Furthermore, the basic 
process in a reflection experiment, multiple elastic scattering, was treated approximately, 
in the so-called Px-approximation, in which the expansion of the elastic scattering cross 
section in spherical harmonics is terminated after the first order [19]. In the medium energy 
range, the higher order phase shifts are by no means negligible [20] and consequently the 
Pi-approximation does not provide a good physical basis to model REELS spectra. 

Nonetheless, this formula is extensively used nowadays and it is commonly believed that 
the "effective" cross section that is extracted from a REELS spectrum in this way indeed 
represents a single scattering loss distribution, although the exact nature of this "effective" 
cross section has never been analyzed. 

Some years later, several authors independently introduced a deconvolution procedure for 
REELS spectra in which elastic scattering is treated more realistically [9, 21]. Although this 
provides a better model to describe REELS spectra, the resulting deconvolution procedure 
leads to results that are very similar to those using Tougaard and Chorkendorffs algorithm. 
This has been attributed to the fact that surface excitations are also neglected in this 
procedure [9, 21, 22]. 

Recently, a deconvolution procedure was introduced that is based on a physical model that 
takes into account the multiple occurence of surface excitations and deflections during elastic 
processes [23]. It is based on the simultaneous deconvolution of a set of two REELS spectra 
taken at different energies or geometrical configurations for which the relative contribution of 
surface and bulk excitations differs substantially. The main assumption in this procedure is 
that the normalized single scattering loss distributions for surface and volume scattering are 
independent of the energy of the probing electron and the angle of surface crossing [7, 23]. 
The essential difference with the earlier algorithms is that the latter are essentially effected 
by reversion of a univariate power series in Fourier space which has no unique solution if 
more than one kind of inelastic proces takes place. The algorithm of Ref.[23] is based on 
reversion of a bivariate power series in Fourier space, which does provide the unique solution 
that is straightforward to interpret quantitatively. 

The procedure based on a bivariate power series reversion was succesfully applied to 
REELS spectra of Si, Cu and Au. The deconvolution procedure is quite involved, how- 
ever, and the question therefore arises how seriously the deficiencies of the earlier methods 
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addressed above affect the results of the deconvolution. The present paper addresses this 
issue. The theory of signal generation in a REELS experiment as well as the different 
deconvolution procedures are reviewed. The reason for the good performance of the Pi- 
approximation for deconvoluting spectra -in spite of the fact that it is not a good model 
to simulate REELS spectra- is elucidated. It is shown theoretically that the procedure of 
Tougaard and Chorkendorff does not result in a single scattering loss distribution. Rather, 
this procedure provides a mixture of surface and volume multiple scattering loss distribu- 
tions and a contribution of mixed terms of any scattering order. Therefore, quantitative 
interpretation of the resulting loss distributions is troublesome. 

Results of analysis of experimental spectra for Fe, Pd and Pt using a bivariate power 
series reversion are compared with results of the earlier algorithms. The bivariate power 
series reversion leads to true single scattering loss distributions that agree quantitatively 
with results based on the theory given by Tung [10] on an absolute scale. The main reason 
for the essential failure of any method based on a univariate power series reversion is that it 
does not eliminate mixed terms corresponding to elecrons that experience both surface and 
volume excitations a certain number of times. 

2. THEORETICAL 

The theory of signal generation in REELS is discussed in detail in Ref [23], here we only 
present a brief synopsis. The flux of medium-energy electrons reflected from a solid surface 
consists of particles that have travelled different pathlengths in the solid. Since the prob- 
ability for (multiple) scattering increases monotonically with the travelled pathlength, the 
outgoing flux consists of groups of electrons that have experienced a different number of in- 
elastic processes. The net energy loss experienced after a given number of inelastic processes 
is subject to fluctuations that can be expressed in terms of a multiple self-convolution of the 
single scattering energy fluctuation distributions as follows from solution of a Boltzmann- 
type kinetic equation [24, 25]. The single scattering fluctuation distributions are given by the 
inverse differential mean free path for scattering. The normalized inverse differential mean 
free path for surface (subscript "s") and bulk (subscript "b") excitations will be denoted by 
w s (T) and Wf,(T) in the following. The normalized bulk mean free path is related to the 
absolute differential mean free path Wb(T, E) via the total (integrated) inelastic mean free 
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path (IMFP) \i(E): 



w b (T,E) = Xi(E)W b (T,E) 



(1) 



The normalized single scattering surface loss distribution is given in terms of the differential 
mean free path integrated over the part of the trajectory in the surface scattering zone, 



where the (dimensionless) quantity (n s (E,6)) represents the average number of surface ex- 
citations in a single surface crossing along the polar direction cos#. The dimensions of 
w s (T, E, 9) and w b (T, E) are the same, reciprocal eV, while the dimensions of the quantities 
W S (T, E, 9) and W b (T, E) are different, as is their physical meaning (see Ref. [23]). It turns 
out [7, 23, 26] that the shape of the single scattering loss distributions depends only very 
weakly on the energy and, in the case of surface excitations, on the surface crossing direction. 
Then the normalized single scattering loss distributions are only a function of the energy 
loss to a good approximation: 



Formulae to establish the differential inverse inelastic mean free path (DIIMFP) Wb(T,E) 
and the differential surface excitation probability (DSEP) W S {T, E, 9) can be found in many 
instances in the literature, in the present work the expressions given in Ref. [10] will be 
employed. The differential mean free path for inelastic scattering is frequently refered to in 
the literature as "cross section" but one should keep in mind that in doing so, the solid state 
character of the medium is completely ignored, which is certainly not justifiable for medium 
energy electrons travelling in solids. 

In the following, the fluctuation distributions after rig-fold surface and n b -fold volume 
scattering will be denoted by wi n ^ and w^, respectively, where the exponent in brackets 
(W) indicate the n — 1-fold selfconvolution of the respective quantity. Multiplying these 
with the relative number of electrons that arrive in the detector after suffering the corre- 
sponding number of energy losses, the so-called partial intensities, a nb>n3 , and summing over 
all scattering orders, one obtains the following expression for the loss spectrum Ul(T): 



W S (T,E,9), as: 



w s (T,E.9) 



W S (T,E,6) 
ME, 9)) 



(2) 



w b (T,E) ~ w b (T) 
w s (T,E,9) ~ w s (T) 



(3) 



oc 



oc 




(4) 



ri(,=0 n s =0 
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where the symbol <g> denotes a convolution over the energy loss T. The spectrum Ul(T) 
represents the loss spectrum divided by the area of the elastic peak and from which the 
elastic peak is removed by choosing a„ ft=0 ns=0 = 0. 

The reduced partial intensities for surface and bulk scattering are uncorrelated: 

^n(,,n s ^rij, X Oi ns (5) 

since the part of the trajectory in the surface scattering zone is rectilinear to a good ap- 
proximation [26]. Although (small) effects of deflections in the surface scattering zone have 
been experimentally observed for a rather pathological case [27], this approach nonetheless 
has proven to constitute an effective approximation [26-29]. 

The reduced bulk partial intensities are given by the well known equation [7, 30] 

oo 

/ W nb {s)Q{s)ds 

a n b = ^ (6) 

JW (s)Q(s)ds 



Here Q(s) is the distribution of pathlengths and W n (s) is the stochastic process for multiple 
scattering, which, in the quasielastic regime, is given by 



/ e \ n p s /^i 

W ^> = (a") — < 7 > 

The partial intensities for surface scattering, A ns , follow Poisson statistics, implying that 
the reduced partial intensities a ns = A n JA ns=Q are given by: 

«» s = — r 8 

with (n s ) = (n l s ) + (n°) representing the sum of the average number of surface excitations 
along the in and outgoing part of the trajectory. Surface excitations along the in- and 
outgoing part of the trajectory are combined in Eqn. (4) by assuming that the shape of 
the norma/izeddistribution of surface energy losses is independent of the surface crossing 
direction 6 [23]. 
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3. DECOMPOSITION OF REELS SPECTRA 



3.1. I: reversion of a univariate power series within the Pi approximation 
(Ref. [6]) . 

In the method of Tougaard and Chorkendorf [6] two approximations are made: (1) surface 
excitations are neglected; and (2) multiple elastic scattering is treated in the so-called Pi- 
approximation. Based on these starting assumptions, a simple deconvolution formula was 
derived that is frequently used also in situations were surface excitations dominate the low 
energy loss region of a spectrum. It is commonly believed that in such cases the method 
provides an "effective" single scattering loss distribution i.e. some kind of linear combination 
of w s (T) and Wb(T) but this has never been proved. 

In the Pi -approximation, the elastic scattering cross section that enters the collision 
term of a Boltzmann-type kinetic equation for the particle transfer, is expanded in spherical 
harmonics and the expansion is terminated after the first order. In this way, a tractable 
analytic solution for the reflection problem has been found [31]. The expression for the 
pathlength distribution Qpi(s) reads: 

Q Pl (s) (xexp(-s/L) (9) 

where the characteristic length L is about twice the transport mean free path L ~ 2\ tr . 
Inserting this result into Equation (6), the bulk partial intensities are found as: 

<\ = (10) 

where k = L/(L + Aj). 

When surface excitations are neglected, Equation (4) simplifies to: 

oo 

y L (T) = J2*nA nb \T) (11) 

n b =0 

Convoluting this equation with a 2 / 'ot\Wb(T) and subtracting the result from Equation (11) 
gives: 

oo 

y L (T) - -w b (T>) ® y L (T - T') = a lWb (T) + V(a n - ^a n _>£ n) (12) 

For a sequence of partial intensities of the form a n = k", all terms in the sum on the right 
hand side vanish. Therefore, within the Px-approximation for elastic scattering, the exact 



solution of Equation (11) is: 

^^(T) = y L (T) - -y^W^T') ® y L (T - T'), (13) 

where Equation (1) has been used. This is the deconvolution formula of Tougaard and 
Chorkendorff for which a simple recursive scheme was proposed by these authors [6]. 



3.2. II: exact reversion of a univariate power series (Ref. [9] and [21]) 

Although Equation (13) is the exact solution of Equation (11) in the Pi-approximation, 
its validity has been questioned by several authors [9, 21], since the Pi -approximation does 
not provide an accurate model for the simulation of REELS spectra. This becomes evident on 
comparing the partial intensities in the Pi-approximation, with more realistic results based 
on the Mott cross section for elastic scattering. Results of such calculations are displayed 
in Figure 1 for electrons travelling in Fe for several energies, along with the corresponding 
Mott cross section [32]. It is seen that the partial intensities in the Px-approximation are 
qualitatively different from the more realistic results based on the Mott cross section. This is 
caused by the influence of the shape of the cross section on the sequence of partial intensities 
[20, 33]. Realizing that the partial intensities represent the number of electrons arriving in 
the detector after a given number of inelastic processes, it is clear that the Px-approximation 
is not a good approach for quantitative simulation of REELS spectra (see Figure 3a below). 

The exact solution of Equation (11) for an arbitrary sequence of partial intensities, not 
necessarily of the form a n = K n , was found independently by Vicanek [9] and the present 
author [21]: 

oo 

w b (T)=J2^yf\T) (14) 

q=l 

The coefficients u q in this solution are functions of the partial intensities, that are derived 
in Appendix I. This procedure was succesfully applied to experimental spectra given by a 
univariate power series in Fourier space, obtained by eliminating multiple volume scattering 
from them prior to application of Eqn. (14) [26, 28]. 

For partial intensities in the Px-approximation (a nb = K nb ), the coefficients u q assume a 
particularly simple form: 

(-1)5+1 

u q = [ -^— (15) 
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Therefore, a non- recursive alternative representation of Eqn. (13) is given by: 

oo 

Kw b (T) = X>l) 9+1 #(n (16) 

q=l 

which is sometimes useful for a theoretical analysis, although Eqn. (13) is simpler for prac- 
tical applications. 

3.3. Ill: exact reversion of a bivariate power series (Ref. [23]) 

The above deconvolution procedures proceed via reversion of a univariate power series in 
Fourier space. However, when more than one type of scattering contributes to the spectrum, 
the spectrum is represented by a bivariate power series in Fourier space and, in consequence, 
a univariate power series reversion does not give the unique solution for the involved mean 
free paths. This problem can be resolved [23] by simultaneously deconvolving a set of two 
spectra taken at different energies or scattering geometries, leading to a different sequence of 
partial intensities. In doing so, it is assumed that the norma/uecfclifferential mean free path 
for volume and surface scattering is independent of the energy and scattering geometry, i.e. 
that the normalized mean free paths in the two spectra are identical [see Eqn. (3)]. This 
turns out to be a very reasonable approximation [7, 23]. Then, the unique solution for the 
mean free path for surface and bulk scattering is given by the expression: 

oo oo 

w b{T) = 

p=0 q=0 

oo oo 

w s(T) = E uS P JA(T') ® yg(T - T') (17) 

where the coefficients u b pq and u s pq are again functions of the partial intensities. The algorithm 
to obtain these coefficients is given in Appendix II. 

4. RESULTS AND DISCUSSION 

In order to compare the procedures described above, they were first applied to results 
of model calculations (as described in Ref. [23]) and subsequently experimental data were 
deconvoluted with the method of Tougaard and Chorkendorff (TC) and with the reversion 
of a bivariate power series. 
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The procedure used to acquire the experimental data has been described in detail before 
[33]. REELS data for Fe, Pd and Pt were measured in the energy range between 300 and 
3400 eV for normal incidence and an off- normal emission angle of 60°, using a hemispherical 
analyzer operated in the constant analyzer energy mode giving a width of the elastic peak 
of 0.7 eV. Count rates in the elastic peaks were kept well below the saturation count rate of 
the channeltrons and a dead time correction was applied to the data. For each material the 
optimum energy combination for the retrieval procedure of two loss spectra was determined 
by inspection of the partial intensities. For most materials the optimum energy combination 
was (1000-3000 eV) [23]. 

All simulated data were calculated for the scattering geometry used in the experiment 
with the method described in Refs. [23, 33]. In a first step, simulated spectra were generated 
neglecting surface excitations and using the bulk partial intensities calculated within the P\- 
approximation. In this case method I exactly returned the model DIIMFP. When the more 
realistic partial intensities based on the Mott cross section are used (see Figure 1), the 
simulated spectra are significantly different. However, the results of method I and II for 
the shape of the retrieved DIIMFP was indistinguishably similar, the only difference is that 
the absolute value of the DIIMFP differs since the first and second order partial intensities 
are different. The reason why both methods give essentially the same results for the shape 
of the DIIMFP is easy to understand on inspection of the partial intensities in Figure 1. 
Although the sequence of partial intensities obtained from MC calculations is qualitatively 
different from the one obtained in the Pi-approximation, the MC-partial intensities for the 
two energies shown are also to a good approximation of the form a n = n n (with a different 
value for k), at least for the first few scattering orders n ~ 1 — 5. Therefore, method I 
essentially results in the correct deconvolution for any spectrum that can be represented 
as a univariate power series of the DIIMFP in Fourier space. This is true for any realistic 
sequence of partial intensities for the first few scattering orders, except for some really 
pathological cases [21]. 

However, when method I is applied to simulated data in which surface excitations are not 
ignored, the results are no longer interpretable in a straightforward way. This can be seen 
in Figure 2 that shows model calculations for loss spectra with plasmon peaks modelled by 
Gaussian peak shapes (15 eV for the surface- and 25 eV for the bulk plasmon). The width 
of the plasmon loss peaks was taken to be 0.4 eV in (a) and 2.4 eV in (b). The retrieved 
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loss distributions clearly contain significant intensity from higher order scattering. 

This becomes understandable by analyzing the loss distribution retrieved by method I. 
In order to do so, the spectrum [Eqn. (4)] is inserted into Eqn. (16), giving: 

kw tc (T) = a lfi w b (T) + a 0A w s (T) - a hl w b {T') <g> w s (T - T') - a 2ja wf\T) + <D(w®) (18) 

where Eqns. (5) and (8) were used. Thus, it is found that the loss distribution of Tougaard 
and Chorkendorff contains a negative contribution equal to the second order mixed and 
surface scattering term. This is clearly observed in Figure 2a. It is also seen that the higher 
orders decrease rapidly with the scattering order. Since the average number of surface exci- 
tations is usually quite small compared to unity, the contribution of surface excitations may 
become negligible for high enough energies, but the essential point is that the second order 
mixed term is always appreciable in magnitude since it remains comparable to the single 
scattering partial intensities for surface (and bulk) scattering. This finding also explains the 
spurious negative excursions for energies corresponding to one surface and one bulk loss in 
the retrieved loss distributions in Ref. [6] and similar works published in the past twenty 
years. 

Similar results were obtained for more realistic simulations for 500 eV electrons in Fe, 
shown in Figure 3. In Figure 3a the dotted curve represents the simulation in the Pi- 
approximation, while the dashed curve is based on Monte Carlo (MC) calculations for the 
pathlength distribution from which the partial intensities are derived via Equation (6). 
The latter are seen to agree better with the experimental data (represented by the noisy 
solid curve). The reason is that the realistic Mott cross section for elastic scattering can 
not be described by a first order Legendre polynomial, as implicitly assumed in the Pi- 
approximation (see Figure la). In Figure 3b, the distribution retrieved with method I is 
shown as the solid curve and is compared with the single surface and volume scattering 
contribution to the spectrum, represented by the chain-dashed and long-dashed curves, 
respectively. The short dashed curve is the total single scattering distribution of the model 
spectrum, i.e. ai^Wb(T) + ao^w s (T). It is seen that the retrieved loss distribution is lower 
than the single scattering loss distribution and that it can not be represented by any linear 
combination of the single scattering contributions. Thus, the method of Tougaard and 
Chorkendorff never gives a smg/escattering loss distribution when surface scattering is not 
negligible. 
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Finally, the method suggested by Chen [22] deserves to be addressed in this context. This 
author proposed to extract the DIIMFP and the DSEP by using the "effective" cross section 
of Tougaard and Chorkendorff obtained from REELS spectra taken at two different energies. 
This method is thus experimentally similar to the bivariate power series reversion, however, 
the results are quite different, as will be shown below. To eliminate the first order bulk 
contribution, giving a retrieved distribution with the first order DSEP as the leading term, 
two Tougaard- Chorkendorff cross sections wf c and w\ c are subtracted from each other 
after one of them is multiplied with v = «i,o//3i,o- Here Oi nb>na are the partial intensities of 
spectrum 1 and f3 nb ,n a are the partial intensities of spectrum 2. This gives: 

k iW I c {T)-k 2 wI c {T) = A lfi w b (T)+A 0A w s (T)-A 1A w b (T')^w s (T-T')-A 2fi w^ (T)+0(w 3 ) 

(19) 

where A nb;7ls = ct n ^ ns — v(3 nb ,n s - It is seen, also in this case, that the second order mixed 
term is always of the same order of magnitude as the first order surface term in the resulting 
(difference) spectrum: 

(20) 

^0,1 

So that in this method the mixed terms also make it impossible to recover the single scat- 
tering loss distributions from two "effective" cross sections obtained from REELS spectra 
taken at two energies (or scattering geometries). The analysis of Chen [22] is quite different 
from the present one. In the work by this author, the mixed terms are "eliminated" by 
replacing the Fourier transform of the bulk DIIMFP (Jj° b (s) in the notation of that paper) 
in Eq.(45) by the same quantity in real space in Eqns. (46)-(51). In this way this procedure 
ignores the mixed terms, but unfortunately, they are of the same order of magnitude as the 
first order surface and bulk term. 

The result of the retrieval using the bivariate power series reversion is compared in Fig- 
ure 3c (dashed lines) with the model DIIMFP and DSEP for 3000 eV (solid lines). In this 
case the retrieved distributions for surface and bulk scattering are separated and are very 
similar to the model distributions. The differences that can be observed are due to the 
assumption that the normalized DIIMFP and DSEP are independent of the energy of the 
probing electron and the direction of surface crossing. Although this is a good approxi- 
mation, it is not exactly true and hence the retrieved distributions agree with the model 
loss distribution (for one energy) within the accuracy of this approximation. If the model 
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DIIMFP and DSEP are taken to be the same for both energies and for the incoming and out- 
going direction of surface crossing, the retrieved distributions exactly agree with the model 
distributions, since Eqn. (17) is the exact solution of Eqn. (4). 

Application of Tougaard and Chorkendorffs method to the 1 keV experimental spectra 
of Fe, Pd and Pt is shown in Figure 4 as open circles. The dashed and dash-dotted curves 
represent the theoretically expected contribution of surface and volume scattering, the solid 
curve is their sum. For all three materials, the experimentally retrieved loss distribution 
in Figure 4 is seen to be significantly lower than the total single scattering contribution 
(solid curves) for energies < 80 eV while for higher energy losses they become similar to 
the theoretically expected result. This makes it clear that the loss distributions retrieved 
by the method of Tougaard and Chorkendorff are not a linear combination of the single 
scattering surface and volume loss distributions being attributable (mainly) to the negative 
second order mixed term, as discussed above. Thus, although the retrieval of these data 
from experimental results is very simple, the usefullness of this procedure is questionable 
since quantitative interpretation of the results is by no means straightforward. 

The method of reversion of a bivariate power series in Fourier space is illustrated in 
Figure 5 for Fe (1000-3000 eV), Pd (500-3400 eV) and Pt (1000-3000 eV). Figure 5a compares 
the measured loss spectra (i.e. the spectra after division by the area of the elastic peak and 
after removal of the latter) as noisy solid curves with simulated spectra (smooth solid curves). 
Reasonable agreement between simulation and experiment is seen in all cases, but significant 
deviations can also be observed. For example, the ionization edges at ~60 eV are not exactly 
reproduced by the simulation in all cases. This is attributable to the fact that it is difficult 
to fit ionization edges with a Drude-Lindhard type expansion of the dielectric function [10]. 

In Figure 5b and c the DSEP and DIIMFP retrieved with the bivariate power series 
reversion are represented by the open circles and compared with result based on the theory 
of Tung and coworkers [10] (solid curves). In all cases the agreement is quite good, except in 
the vicinity of ionization edges. Furthemore, the main feature in the DIIMFP of Fe is a bit 
sharper than predicted by theory and the DSEP for Pd is significantly higher in the energy 
range between 10-30 eV. Note that the retrieved distributions were not scaled in any way, 
the normalized loss distributions are compared on an absolute scale. 

The DSEP always exhibits a negative excursion for an energy corresponding approxi- 
mately to the average volume loss. This feature, that has a clear physical explanation [26] , 
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is also seen in the theoretical results and is a consequence of the coupling of surface and 
bulk excitations that are essentially different (orthogonal) modes of the inelastic interaction. 
The coupling term is included in the DSEP, as defined in the present work [23, 27], and it 
is always negative. So the DSEP is in fact a difference of two probabilities, the pure surface 
excitation term and the coupling term that takes into account the orthogonality of the sur- 
face and bulk modes, and the resulting DSEP exhibits a negative excursion whenever the 
coupling term exceeds the pure surface term. 

Finally, in Figure 6 the present results for the volume single scattering loss distribution 
is compared with results based on optical data from various sources in the literature [1, 2, 
10, 34]. The literature data exhibit significant discrepancies. For Pd, the present results 
agree well with those of Tung [10] for energies below ~50 eV. For Pt, the present results 
excellently match the data by Tanuma [34] and Palik [1]. In both cases, the ionization edge 
is more pronounced in the present data. This feature is well reproduced by first principles 
calculations based on density functional theory calculations beyond the ground state for 
both materials (dash-dotted lines) [2]. 



5. SUMMARY AND CONCLUSIONS 



The single scattering loss distribution for inelastic electron scattering was extracted from 
experimental data for Fe, Pd and Pt. This was done by means of a recently proposed proce- 
dure [23] in which it is assumed that the normalized differential mean free path for inelastic 
scattering (DIIMFP) does not depend on the energy of the probing electron. Furthermore, 
it is assumed that the normalized differential surface excitation probability (DSEP) does not 
depend on the surface crossing direction and the primary energy. Under these assumptions, 
a set of two spectra with a different sequence of partial intensities can be simultaneously 
deconvoluted to give the normalized true single scattering loss distributions. The decon- 
volution procedure proceeds via the reversion of a bivariate power series in Fourier space. 
Experimentally, two spectra with a different sequence of partial intensities can be obtained 
by measuring the loss spectrum at two different energies, as studied in the present work, or, 
alternatively, for two different scattering geometries. The latter case, which was not studied 
here, is a bit more demanding experimentally, but has the great advantage that the DIIMFP 
and DSEP of the two spectra are indeed identical. In other words, it is then only the angu- 
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lar dependence of the DSEP (for the in- and outgoing beam) that may slightly distort the 
deconvolution results. However, the present results based on spectra taken at two different 
energies also compare very well with theoretical results since the energy dependence of the 
normalized single scattering loss distributions is quite weak in reality. 

A detailed comparison was made with earlier methods based on reversion of a univari- 
ate power series. The most well known and frequently used deconvolution procedure by 
Tougaard and Chorkendorff [6] belongs to this class of procedures. It was shown, both the- 
oretically and on the basis of model calculations, that when such a method is applied to a 
spectrum described by a bivariate power series (in the surface and volume single scattering 
loss distribution), it results in a not very well defined mixture of multiple volume, surface 
and mixed scattering terms. In particular, the second order mixed term (corresponding to 
electrons that experience one surface and one volume energy loss) gives a strong negative 
contribution in the resulting loss distribution. This explains the spurious negative excursion 
at energies corresponding to one bulk and one surface loss in spectra deconvoluted by the 
method of Tougaard and Chorkendorff that have been published over the past twenty years 
(see e.g., Ref. [6]). 

Finally, and most importantly, the above findings confirm the theoretical model for 
REELS based on a bivariate power series in the DSEP and the DIIMFP in Fourier space 
[23] with unprecedented detail. 
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7. APPENDIX 



7.1. Reversion of a univariate powers series 

We consider the loss spectrum y L written as a power series of the variable w in Fourier 
space (i.e. a multiple self-convolution in real space): 

oo 

y L (T) = $>V n) (T) (21) 

n=0 

where a = since it is assumed that the elastic peak has been removed from the spectrum. 
Since the leading term in y L is the first power of w, the solution must be of the form: 

oo 

w{T) =Y. u MT) iq) ■ (22) 

q=l 

Inserting Eq.21 into this ansatz one finds: 

oo oo , < 

M^) = E^{^ a ^ (n)(T) }' ? ' (23) 

q=l n=l 

and, equating coefficients of equal powers in the above equation: 

1 = 

= M 1 a 2 + u 2 aicei 

= + M2(q;i«2 + ct20ti) + uzaiaiai 

= UiCti + U2{ot20t2 + «itt3 + Q-zCti) + u^{a20c\CXi + aiQ^i + o.\aia 2 ) + 1x40:10:10:10:1 



the solution for the unknown coefficients u q is easily seen and can be found in any standard 
textbook [35, 36] in a slightly different and more explicit albeit less general form. The 
interesting point in the present notation, however, is that the subscripts in the coefficients 
of u q in the above equations represent the sum of all permutations of the partitions of n (the 
order) with q factors in them. 
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7.2. Reversion of a bivariate power series. 

When two loss spectra with a different sequence of partial intensities are measured, 

oo oo 

n s =0 n;,=0 
oo oo 

VlAT) = £ £ ® w^\T - T'), (24) 

n s =0 rtj,=0 

with a ,o = A),o = 0, a unique solution for the unknowns w s (T) and Wb(T) can be established 
by reversion of this bivariate power series (in Fourier space). 

Formally, the reversion of this bivariate power series is effected by the expansion: 



oo oo 



p=0 q=0 

oo oo 

w s(T) = «rfiCn ® VlI(T - T') (25) 

p=0 q=0 

with Uq = Mq = 0. This can be seen by substituting Equation (24) back into Equation (25) 
and equating coefficients of equal powers of Wb and w s . This gives the equations for the 
unknown coefficients u h p and u s p . 

In doing so, one is faced with the problem of evaluating the convolution of y^iT) with 
Vl\(T) that can be expressed as: 



oo oo 



vigCO ® yg(T -?*) = ££ 7 M ,n 6 ,n.«;S B6) (T') ® ^>(T - T') (26) 

n s =0 «6=0 

The coefficients r y Ptq ,n b ,n s are found to be given by: 

lp,q,n b ,n s = (dp, Pq)^^) (27) 

These coefficients are equal to the sum of all possible terms with p factors in etk,i and q factors 
in f3 m>n , whose indices "add up" to the target index combination (n s ,rif,) = (k + m,l + n). 
For example: 

( a V(2,l) = a l,0«l,l + tt0.1«2,0 + Q!l,lQ!l,0 + tt2,0«0,l 

(ai,/3 2 )(!_2) = «o,i^o,iA,o + tto,iA,o/3o,i + «i,o/5o,i/5o,i (28) 

Obviously, one has 

Tiw.rn.n, =0 for all p + q > n s + n b (29) 
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since the target index combination (n s , rib) can om y be produced by a number of factors less 
than or equal to (n s + rib) when a ,o = A),o = [see Equation (24)]. Furthermore, one has: 

7l,0,n b ,n s = a n b ,n a 

= Pn b ,n s (30) 

Based on these guidelines, the coefficients r y Pt q t n b ,n s can readily be established by means of a 
recursive algorithm for any sequence of partial intensities Oi nb>na and f3 nb ,n s - 

Inserting Equation (24) into Equation (25) and equating coefficients, one finds a set of 
two equations with two unknowns for the first order bulk coefficients: 



OO CO 



p=0 q=0 

OO OO 

= Yl Yl = «o,iA),i + u b lfi ao,i (31) 

p=0 q=0 



The solution is: 



<o = 



/3o,iai j0 — /3i,o a o,i 



«0,lPl,0 — tti,oPo,i 



Similarly, for the first order surface coefficients one finds: 



M l,0 



/3i,o«o,i — Auc*i,o 



"« = „ a ""*„ a < 33) 



The equation determining the higher order {p + q > 1) coefficients 



OO OO 



^ ] ^ ] u p,qlp,q,n b ,n 3 (34) 
p=0 g=0 

can be split into two parts by using property (29), changing the order of the summation and 
writing: 

n s +n b — 1 p n s +n b 

= Uq,p-qlq,p-q,n b ,n s + ^p,n s +n;,-p7p,n 6 +n s -p,n;„n s (35) 

p=0 g=0 p=0 

The first term represents a (n s + + l)-dimensional vector containing the coefficients 
u p<n b>q <n s i which have been established during the previous step of the algorithm. The latter 
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term is the product of a square matrix with the same dimension and the unknown vector 
(u ,n b +n s , ui,n 6 +n.-i, M 2 ,n 6 +n s -2, ...u nb+nsfi ). Thus, for each value of n b + n s , the corresponding 
coefficients are obtained by solution of a system of n b + n s + 1 linear equations with as many 
unknowns. Consecutively performing this procedure for values of the total scattering order 
n b + n s = 1, 2, ....n max , where n max is the collision order where convergence of the series 
Equation (25) is attained, leads to the desired coefficient matrices u b p q and u s p q . Note that 
Eqns. (34) and (35) are identical for the surface and bulk coefficients, only the calculation 
for the first order term is different [see Eqns. (32) and (33)]. Although faster algorithms 
for reversion of multivariate power series are available [37-39], these are quite complex and 
were not considered here since convergence of the proposed algorithm is typically attained 
for n max < 7 for energy loss ranges extending up to 100 eV. 
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FIG. 1: (a) Differential Mott cross section for elastic scattering for Fe for two electron energies, (b) 
Reduced partial intensities a nf> for volume scattering for several energies in Fe. Solid curves: Monte 
Carlo calculations based on the Mott scattering cross section. Dashed curves: Pi-approximation 
(Equation (10)) 
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FIG. 2: Simulated spectra (solid lines) for a model DIIMFP and DSEP represented by a Gaussian 
peak shape with a mean energy loss of 25 and 15 eV, respectively. The dashed curves are the result 
of application of Tougaard and Chorkendorff's algorithm Eqn. (13). It is seen that the latter does 
not eliminate the higher order terms, (a.) a model plasmon width of 0.4 eV was used, (b) a model 
plasmon width of 2.4 eV was used. 
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FIG. 3: (a) experimental REELS spectrum for 3000 eV-electrons backscattered from an Fe-surface 
(solid curve), spectrum simulated with Equation (4) based on the Mott cross section for elastic 
scattering (dashed curve) and simulated spectrum within the Pi-approximation for elastic scatter- 
ing (dotted curve), (b) loss distribution retrieved from the simulated spectrum in (a) based on the 
Mott cross section for elastic scattering using the Method of Tougaard and Chorkendorff [6] (solid 
curve) . The long-dashed curve represents the single volume scattering contribution of the simulated 
spectrum, the dash-dotted curve is the contribution of single surface scattering. The short-dashed 
curve is the total single scattering contribution in the spectrum in (a), (c.) normalized DIIMFP 
and DSEP retrieved from the simulated spectra in Figure 5(a) using the bivariate power series 
reversion [23] (dashed curves) compared with the DIIMFP and DSEP used in the simulation (solid 
curves) . 
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FIG. 4: Loss distribution obtained with the method of Tougaard and Chorkendorff [6] (open circles) . 
The dash-dotted and dashed curves are the normalized DIIMFP and DSEP calculated with the 
formulae in [10] and multipled with the according single scattering partial intensities, the solid 
curve represents their sum. In all cases, it can be seen that the retrieved spectrum is nota. linear 
combination of the single scattering terms, (a) Fe; (b) Pd; (c) Pt. 
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FIG. 5: (a) Measured energy loss spectra for Fe, Pd and Pt for the energies indicated in the 
legend after removal of the elastic peak (noisy curves). The smooth curves are the simulated loss 
spectra using theoretical shapes for the DSEP and DIIMFP. (b) Retrieved normalized DSEP (open 
circles), (c) Retrieved normalized DIIMFP (open circles). The solid curves in (b) and (c) represent 
the result based on the theory by Tung and coworkers [10]. 
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FIG. 6: Comparison of the DIIMFP obtained in the present work (open circles) with theoretical 
calculations using optical data taken from different sources [1, 2, 10, 34] as indicated in the legend, 
(a) Pd, 1000 eV. (b) Pt, 1000 eV. 
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